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\ Abstract 

In this paper, new travelling wave solutions to the Fornberg-Whitham equation 

i Ut U X xt U x ~\~ UU X — UU X xx ~\~ Stt^Ujj; 

d 

are investigated. They are characterized by two parameters. The expresssions of the 
periodic and solitary wave solutions are obtained. 
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Recently, Ivanov [1] investigated the integrability of a class of nonlinear 
^ \ dispersive wave equations 

ut - u xxt + 8 x (ku + au 2 + (3u z ) = vu x ii xx + iuu xxx , (1.1) 



where and a, (3, 7, k, v are real constants. 
The important cases of Eq.(l.l) are: 
The hyperelastic-rod wave equation 
u t - Uxxt + 3uu x = ^(2u x u xx + uu xxx ), (1.2) 



* Corresponding author. 

Email address: zhoujiangbo@yahoo.cn (Jiangbo Zhou). 

Preprint submitted to 6 August 2009 



has been recently studied as a model, describing nonlinear dispersive waves 
in cylindrical compressible hyperelastic rods [2]- [7]. The physical parameters 
of various compressible materials put 7 in the range from -29.4760 to 3.4174 
[2,4]- 

The Camassa-Holm equation 

u t ~ u xxt + 3uu x = 2u x u xx + uu xxx , (1-3) 



describes the unidirectional propagation of shallow water waves over a flat 
bottom [8,9]. It is completely integrable [1] and admits, in addition to smooth 
waves, a multitude of travelling wave solutions with singularities: peakons, 
cuspons, stumpons and composite waves [9]-[12]. The solitary waves of Eq.(1.2) 
are smooth if n > and peaked if k — [9,10]. Its solitary waves are stable 
solitons [13,14], retaining their shape and form after interactions [15]. It models 
wave breaking [16]- [18]. 

The Degasperis-Procesi equation 

u t - u xxt + 4uu x = 3u x u xx + uu xxx , (1.4) 



models nonlinear shallow water dynamics. It is completely integrable [1] and 
has a variety of travelling wave solutions including solitary wave solutions, 
peakon solutions and shock waves solutions [20]-[27]. 

The Fornberg-Whitham equation 

Ut t^xxt ~t~ U x ~\~ uu x = uu xxx -\- <5u x u xx , (1-5) 



appeared in the study qualitative behaviors of wave-breaking [28]. It admits 
a wave of greatest height, as a peaked limiting form of the travelling wave 



solution [29], u(x,t) = v4exp(— | 
It is not completely integrable [1] 



x 3 t 



where A is an arbitrary constant. 



The regularized long-wave or BBM equation 

u t - u xxt + u x + uu x = 0, 



;i.6) 



and the modified BBM equation 



Ut - u xxt + u x + 3u u x = 0, 



;i-7) 



have also been investigated by many authors [30]- [38]. 
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Many efforts have been devoted to study Eq.(1.2)-(1.4),(1.6) and (1.7), 
however, little attention was paid to study Eq.(1.5). In [39], we constructed two 
types of bounded travelling wave solutions «(£)(£ = x — ct) to Eq.(1.5), which 
are defined on semifinal bounded domains and called kink-like and antikink- 
like wave solutions. In this paper, we continue to study the travelling wave 
solutions to Eq.(1.5). Following Vakhnenko and Parkes's strategy in [19], we 
obtain some periodic and solitary wave solutions w(£) to Eq.(1.5) which are 
defined on (— oo,+oo). The travelling wave solutions obtained in this paper 
are obviously different from those obtained in our previous work [39]. To the 
best of our knowledge, these solutions are new for Eq.(1.5). Our work may help 
people to know deeply the described physical process and possible applications 
of the Fornb erg- Whit ham equation. 

The remainder of the paper is organized as follows. In Section 2, for com- 
pleteness and readability, we repeat Appendix A in [19], which discuss the 
solutions to a first-order ordinary differential equaion. In Section 3, we show 
that, for travelling wave solutions, Eq.(1.5) may be reduced to a first-order 
ordinary differential equation involving two arbitrary integration constants a 
and b. We show that there are four distinct periodic solutions corresponding 
to four different ranges of values of a and restricted ranges of values of b. A 
short conclusion is given in Section 4. 



2 Solutions to a first-order ordinary differential equaion 

This section is due to Vakhnenko and Parkes (see Appendix A in [19]). For 
completeness and readability, we repeat it in the following. 

Consider solutions to the following ordinary differential equation 



and tpi, (f2, <fi3, fi are chosen to be real constants with ip\ < if 2 < V 9 ^ ^3 ^ 
if 4 . 

Following [40] we introduce ( defined by 



(m) 2 = e 2 /M, 



(2.1) 



where 



= {<P- ¥>i)(¥> - ^2)^3 - <f)(<f4 ~ <f), 



(2.2) 



(2.3) 
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so that Eq.(2.1) becomes 
(^C) 2 = /(¥>)• 



(2.4) 



Eq.(2.4) has two possible forms of solution. The first form is found using 
result 254.00 in [41]. Its parametric form is 



_ (p 2 — (pinsn 2 (w\m) 
^ 1 — nsn 2 (w\m) 



£ = ^(wif! + (<p2 - ¥>i)n(n; w\m)), 



(2.5) 



with w as the parameter, where 



m= 7 w t,P= ttV^ 4- ^ 2 )(^3-^i),w = K, (2-6) 



and 



n 



V?3 - V2 



(2.7) 



In (2.5) sn(w\m) is a Jacobian elliptic function, where the notation is as used 
in Chapter 16 of [42]. U(n;w\m) is the elliptic integral of the third kind and 
the notation is as used in Section 17.2.15 of [42]. 

The solution to (2.1) is given in parametric form by (2.5) with w as the 
parameter. With respect to w, <p in (2.5) is periodic with period 2K(m), where 
K(m) is the complete elliptic integral of the first kind. It follows from (2.5) 
that the wavelength A of the solution to (2.1) is 

A = — 1^(70) + (<p 2 -<pi)n(n\m)\, (2.8) 
ep 



where U(n\m) is the complete elliptic integral of the third kind. 

When </? 3 = if 4, m = l, (2.5) becomes 

_ (f2 — fin tanh_ w 
^ 1 — ntanh w ' (2.9) 
£ = I(» _ 2tanh~ 1 ( v ^tanhw)). 



The second form of solution of (2.5) is found using result 255.00 in [41]. Its 
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parametric form is 



_ ip 3 — (^ 4 nsn 2 (w|m) 
^ 1 — nsn' 2 (w\m) ' ^2 10) 

£ = ^p( WL Pi ~ (^4 - ^3)n(n; w\m)), 



where m,p,w are as in (2.6), and 

n= ^~^ . (2.11) 
tp 4 - ip 2 



The solution to (2.1) is given in parametric form by (2.10) with w as the 
parameter. The wavelength A of the solution to (2.1) is 

2 

A = — \(piK(m) - (ip 4 - <p 3 )Il(n\m)\. (2.12) 
sp 



When if i = tpz, m — 1, (2.10) becomes 

_ <£s — ip^n tanh 2 w 

1 — ntanh w ' (2.13) 
f = 1(^2 +2 tanh- 1 ( v ^ tanh w)). 



3 Periodic and solitary wave solutions to Eq.(1.5) 

Eq.(1.5) can also be written in the form 
[u t + uu x ) xx = u t + uu x + u x . (3.1) 



Let u = if(^) + c with £ = x — ct be a travelling wave solution to Eq.(3.1). 
then it follows that 

= m + ^- ( 3 - 2 ) 



Integrating (3.2) twice with respect to £, we have 

(m) 2 = ^ 4 + !^ 3 + V + &), (3.3) 



where a and 6 are two arbitrary integration constants. 
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Eq.(3.3) is in the form of Eq.(2.1) with e = \ and f(<p) = (ip* + §y? 3 + 



atp 2 + b). For convenience we define g(<p) and h(ip) by 

g 

f(f) = ¥> 2 0(¥>) + &, where 9(f) =f 2 + ^f + a, (3-4) 

/'(<p) = 2(ph((p), where /i(^) = 2ip 2 + Atp + a, (3.5) 

and define <^l, <^r, &£, and Or by 

^L = -^(2 + v / 4^),^ = -^(2- v / 4^), (3.6) 

6 L = -<plg(tp L ) = j - 2a + ~ + ?(2 - a)^!^, (3.7) 

^ = -<p%g(tp R ) = j - 2a + I - H( 2 - a ) v / 4^. (3.8) 



Obviously, </?£,, are the roots of h(<p) = 0. 

In the following, suppose that a < 2 and a ^ such that f(f) has three 
distinct stationary points: ipi, fa, and comprise two minimums separated by 
a maximum. Under this assumption, Eq.(3.3) has periodic and solitary wave 
solutions that have different analytical forms depending on the values of a and 
b as follows: 

(1) a < 

In this case <fL < < <pr and f(fi) < fi.fR)- For each value a < 
and < b < &r (a corresponding curve of f(ip) is shown in Fig. 1(a)), there 
are periodic inverted loop-like solutions to Eq.(3.3) given by (2.5) so that 
< m < 1, and with wavelength given by (2.8); see Fig. 2(a) for an example. 

The case a < and b = 6r (a corresponding curve of f(f) is shown in 
Fig. 1(b)) corresponds to the limit p^ = fi = fR so that m — 1, and then the 
solution is an inverted loop-like solitary wave given by (2.9) with ip 2 < <p < <Pr 
and 

<P! = --(2 + 3V4 - 2a + 2\/4 + Gv 7 !^^), (3.9) 
6 

^ 2 = -1(2 + 3V4 -2a- 2y[l + Gv 7 !^^), (3.10) 
6 

see Fig. 3 (a) for an example. 
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Fig. 1. The curve of f(<p). (a) a = -50, b = 233; (b) a = 
a = 1.5, 6 = -0.05; (d) a = 1.5, b = 0; (e) a = ^, 6 = -0.1; (f) a = ^ 
a = 1.9, 6 = -0.24; (h) a = 1.9, b = -0.2010. 

(2) < a < f 



50, b = 374.1346; (c) 
= 0; (g) 



-1 n 

In this case tpL < < and /(<^l) < /(0). For each value < a < -gp 
and bft < b < (a corresponding curve of is shown in Fig. 1(c)), there 
are periodic hump-like solutions to Eq.(3.3) given by (2.5) so that < m < 1, 
and with wavelength given by (2.8); see Fig. 2(b) for an example. 

The case < a < and b = (a corresponding curve of f(f) is shown in 
Fig. 1(d)) corresponds to the limit y?3 = <y9 4 = so that m = 1, and then the 
solution can be given by (2.9) with ip x and ip 2 given by the roots of g((p) = 0, 
namely 



1 



4 1 



tp 1 = -~- -V16 - 9a, y? 2 = -- + -a/16 - 9a. 



(3.11) 



In this case we obtain a weak solution, namely the periodic upward-cusp wave 
<p = ^ - 2j£ m ), (2j - l)£ m < £ < (2j + l)£ m , j = 0, ±1, ±2, ■ ■ • , (3.12) 



where 

<p(0 = (^ 2 -^itanh 2 (e/4))cosh 2 (e/4), (3.13) 

and 

e m = 4tanh" 1 (3.14) 

V <pi 
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see Fig.3(b) for an example. 
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Fig. 2. Periodic solutions to Eq.(3.3) with < m < 1. (a) a = —50, b 
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Fig. 3. Solutions to Eq.(3.3) with m = 1. (a) a = -50, b = 374.1346; (b) a = 1.5, 
b = 0; (c) a = ^, b = 0; (d) a = 1.9, b 

(3) a 
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1 6 

In this case tpi < <pn < and /(<£l) = /(0). For a = -g^ and each value 
6r < 6 < (a corresponding curve of /(</?) is shown in Fig. 1(e)), there are 
periodic hump-like solutions to Eq.(3.3) given by (2.10) so that < m < 1, 
and with wavelength given by (2.12); see Fig. 2(c) for an example. 
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The case a = and b = (a corresponding curve of /(</?) is shown in 

Fig. 1(f)) corresponds to the limit v 9 i = V 9 2 = V 9 l = — § an d y?3 = y?4 = so 
that m = 1. In this case neither (2.5) nor (2.10) is appropriate. Instead we 
consider Eq.(3.3) with /(</?) = ^y? 2 ((^ + ^) 2 and note that the bound solution 

has — < <£> < 0. On integrating Eq.(3.3) and setting (f = at £ = we 
obtain a weak solution 

^ = ^xp(-i^|)-| (3.15) 

i.e. a single peakon solution with amplitude ^, see Fig.3(c). 

(4) f < a < 2 

In this case tpL < <Pr < and /(</?l) > /(0). For each value 4^ < a < 2 and 
b R < b < b L (a corresponding curve of /(</?) is shown in Fig. 1(g)), there are 
periodic hump-like solutions to Eq.(3.3) given by (2.10) so that < m < 1, 
and with wavelength given by (2.12); see Fig. 2(d) for an example. 

The case 4jp < a < 2 and b = b L (a corresponding curve of f((f) is shown 
in Fig. 1(h)) corresponds to the limit (fii = <y9 2 = (fi so that m = 1, and then 
the solution is a hump-like solitary wave given by (2.13) with <p^ < <p < ^3 
and 

y> 3 = -(-2 + 3V4 - 2a - 2^4 - 6 V / 4^ = 2T), (3.16) 
6 

^4 = -(-2 + 3V4 - 2a + 2\/4 - 6 v / 4^a), (3.17) 
6 

see Fig.3(d) for an example. 

On the above, we have obtained expressions of parametric form for periodic 
and solitary wave solutions </?(£) to Eq.(3.3). So in terms of u — <p(£)+c, we can 
get expressions for the periodic and solitary wave solutions u(£) to Eq.(1.5). 



4 Conclusion 

In this paper, we have found expressions for new travelling wave solutions 
to the Fornberg-Whitham equation. These solutions depend, in effect, on two 
parameters a and m. For m — 1, there are inverted loop-like (a < 0), single 
peaked (a = and hump-like < a < 2) solitary wave solutions. For 
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m = l,0<a<-gpOrO<m<l,a<2 and a ^ 0, there are periodic 
hump-like wave solutions. 
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